Trapping and cooling of rf-dressed atoms in a quadrupole magnetic field 
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We observe the spontaneous evaporation of atoms confined in a bubble-like rf-dressed trap [l|]. 
The atoms are confined in a quadrupole magnetic trap and are dressed by a linearly polarized rf 
field. The evaporation is related to the presence of holes in the trap, at the positions where the rf 
coupling vanishes, due to its vectorial character. The final temperature results from a competition 
between residual heating and evaporation efficiency, which is controlled via the height of the holes 
with respect to the bottom of the trap. The experimental data are modeled by a Monte-Carlo 
simulation predicting a small increase in phase space density limited by the heating rate. This 
increase was within the phase space density determination uncertainty of the experiment. 



PACS numbers: 32.80.Lg, 03.75.-b, 32.80.Pj 
I. INTRODUCTION 



Ultracold atoms in versatile traps is a subject of 
extensive studies. Developing alternative trap geome- 
tries enables the exploration of new physical situations. 
Amongst the most remarkable results obtained with 
atoms in non harmonic traps, one may cite the observa- 
tion of superfluid to Mott insulator transition in optical 
lattices Josephson oscillations in a double well 
or trapping in a ring geometry Q. Recently, radio- 
frequency (rf ) fields were used together with a static mag- 
netic field to produce a quasi 2D trap [J H, Q and a dou- 
ble well 01 • Other trapping geometries were proposed, 
based on this promising rf-dressing technique, such as 
rings d, [§| or lattices [l^l- One presents here the im- 
plementation of rf-dressing inside a quadrupole magnetic 
trap. 



Quadrupole traps are usually not popular for ultra- 
cold atom trapping due to Majorana losses at the centre, 
where the field vanishes. The losses are avoided in TOP 
traps by adding a rotating uniform magnetic field, at a 
frequency in the 10 kHz range flT'l. In the trap presented 
in this paper, the magnetic field oscillates at a few MHz, 
resulting in an rf-dressed quadrupole trap. The atoms 
are then located away from the region of zero magnetic 
field. Still, the rf polarization is responsible for the ex- 
istence of regions of zero rf coupling. This results in a 
leakage of the higher kinetic atomic population. In this 
Letter, we give evidence for a temperature decrease and a 
non exponential atom number loss, characteristic for an 
evaporation process, for different experimental parame- 
ters. 
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II. TRAP DESCRIPTION 

The trap relies on rf coupling between Zeeman sub- 
levels in an inhomogeneous static magnetic field. The 
basic idea of the rf-dressed potentials was first proposed 
by Zobay and Garraway and experimentally demon- 
strated a few years later |5I| . We will recall here the main 
features of this trap and develop further the consequences 
of an inhomogeneous rf coupling. 

A static trapping magnetic field i?(r)eB(r) present- 
ing a local minimum Bq is used together with a rf os- 
cillating magnetic field i?rf cos(a;rfi) %■ The rf field is 
detuned by A = Wrf — from the coupling frequency 
= gLfJ'sBo/h at the static magnetic field minimum, 
where (Jl is the Lande factor and fiB the Bohr magne- 
ton. The linearly polarized rf field couples the Zeeman 
substates \F,mp) and \F, m'p) with m'p = to^^ ± 1 at the 
positions where the magnetic splitting is close to the rf 
frequency. In the limit of a large coupling, this results in 
a dressing of the mp levels into adiabatic states whose 
energy are represented on fig. [TJ In the following, we will 
concentrate on the upper adiabatic state, undergoing a 
potential 

Vd{v) = F {{Vb{v) - huj^tf + h^n\v))"^ . (1) 

This expression is obtained within the rotating wave ap- 
proximation. Here, mp — F ~ 2m the case of ®^Rb in its 
5si/2 ground state. VB(r) = gi^B-B(r) is the magnetic 
energy shift between two adjacent Zeeman sublevels, and 

f}(r) = 17isin(ey,eB(r)) (2) 

is the Rabi frequency of rf coupling, of maximal ampli- 
tude ill = gLfJ'BBri/{2h,). ^ll is independent of r to a 
very good approximation. As deduced from this formula, 
the rf coupling vanishes at the points where the static 
field is parallel to the radio frequency polarization Gy, 
that is along the y axis. We will see later on that in the 
case of a quadrupole static magnetic field, these points 
appear as two "holes" on the equator of our bubble-like 
trap. 
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FIG. 1: Energy of the dressed levels of Rb in the magnetic 
quadrupole trap described in the paper, plotted along the 
vertical coordinate in the vicinity of the potential minimum 
at z — —ro- The parameters are Wrf = 2it x 3.1 MHz, b' — 
150 G/cm and i}i — 2tt x 100 kHz. The five dressed sublevels 
for a F = 2 spin state are plotted, as well as two bare states 
for comparison, mp = —1, dashed, and mp = 2, dash-dotted. 
The atoms are trapped in the upper dressed potential mp ~ 
F — 2. Gravity was taken into account. 



The simplest situation occurs when the rf coupling 
ri(r) is uniform. Then, for a positive detuning A > 0, 
the dressed potential Vd(r) is minimum all over an iso- 
B surface, defined by Vb(i") = gLfJ-BB{r) — hiVrf- This 
is no longer true if fl varies in space, due to the spa- 
tial dependence of the orientation of the static field: this 
breaks the invariance on the iso-B surface, leading to a 
finite number of minima related to the minima of r2(r). 
On the other hand, gravity tends to push the atoms at 
the bottom of the iso-J5 surface. To describe correctly 
the potential geometry, we thus need to take both effects 
into account. As we will see, the potential minima will 
not always sit exactly on the iso-_B surface, but they will 
be close to it in most relevant experimental cases. 

We concentrate in the following on the case of a 
quadrupole trap with an horizontal symmetry axis. This 
axis is labeled y in our experimental setup, see fig. [21 
and is the same as the rf polarization. Let b' be the mag- 
netic field gradient in the radial direction of coordinate 
p = {x^ + z^y^'^, and let us define a as a = gLfJ-sb'/h. 
One has B(r) = b'{pep — 2yey) and therefore ri(r) = 
ilip/ \/ p"^ + 4y2. In this situation, an atom in the upper- 



most dressed state mp = F undergoes the potential 
Vtot {x, y, z) = mgz 



(3) 



-I- Fh^ {a^Jp^ + 4y2 - uo^if + 



p- 



p^ -f Ay 



where m is the atomic mass and mgz the gravitational 
potential. The relevant iso-i? surface is then an ellip- 
soid of equation p^ + Ay^ = Tq, where rp is related to 
through ro — uj^-f/a. The rf coupling vanishes along 
the y axis, where the magnetic field is parallel to the rf 
polarization. 

To visualize the potential, let us first remark that the 



Quad 2, 




FIG. 2: Simplified scheme of the trap. After accomplishing 
evaporative cooling in the QUIC trap (Quad -I- loffe coils) 
centered at position A, the atoms are loaded into the dressed 
rf trap (Quad + loffe + rf coils). The current in the loffe coil 
is then slowly ramped down to zero for transferring the atoms 
into the dressed trap based on the quadrupole field of axis y 
(Quad -I- rf coils). The rf field polarization is aligned along 
the axis of the quadrupole. 



extrema are located in the vertical plane x — 0. In this 
plane, p = \z\. One can rewrite Vtot(j/, z) in units of mgrQ 
as: 



Vtot (a; = O,y,z)/{mgro) 



(4) 



ro \ 



y/ z'^ + Ay"^ 
ro 



^ 



The parameter f3 — FTia/mg is the ratio of the magnetic 
and the gravitational force. It should be larger than 1 
for the magnetic trap to compensate gravity and in our 
experiment (3 ~ 10. Therefore, the behaviour of the sec- 
ond term in the potential expression is dominant. Indeed, 
the exact study of the minima position of Vtot shows than 
for /3 ^ 1, the minima lie very close to the iso-B surface 



p^ + Ay^ 



To discuss qualitatively the potential ge- 



ometry, we will thus assume that the atoms sit exactly 
on this ellipsoid, which cancels the first term under the 
square root. We are left with: 

z \z\ f2i 

Vtot{z) / {mgro) ~ h /3 , with \z\ < ro . (5) 

ro ro ujii 

Depending on the ratio lu,-[/D,i of the rf frequency to 
the rf coupling, two different situations may occur: for 
cjrf/^^i < /3, or equivalently mgro < Fhfli, the coupling 
inhomogeneity is the dominant effect, and the potential 
minima sit at the points where the coupling term is the 
smallest. This corresponds to the positions where p = 0, 
that is at the intersections of the iso-B ellipsoid and the 
y axis (see fig.IHj). Note that at these points, the rf cou- 
pling vanishes, leading to spin flip losses into untrapped 
states. In the following, we will refer to these two points 
as to the "holes". On the other hand, for ujj-i/Ui > /3, or 
mgrQ > Fhfli, gravity dominates and there is a single 
potential minimum at the very bottom of the ellipsoid, 
far from the holes (fig. In this single well, and in 
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FIG. 3: Contour plot of the potential for Rb atoms in the yz plane, for different values of the rf coupling relative to uj-d. 
The parameters are uj^{/2it = 3.1 MHz and a magnetic gradient of 150 G/cm in the z direction, such that (3 — 9.85 and 
uiif/P = 2-K X 305 kHz. (a) large coupling uJif/Qi < 13 {Q,i/2n — 1 MHz); the two minima are located at the position of the 
holes, (b) intermediate coupling Wrf/fii = /3 (f2i/27r = 305 kHz), (c) small coupling cjrf/fii > l3 {Q,i/2tv = 100 kHz); a single 
minimum is present, at the bottom of the shell. 



the harmonic approximation, the oscillation frequencies 
up to first order in are: 

= \f^- (6) 
V '^o 

ujy = 2J^Jl i = 2uj,Jl~ (3^ . (7) 

V ro y mgro V uj^f 

The vertical oscillation frequency was derived in [l| . One 
easily identifies the pendulum frequency along x. Along 
y, the pendulum frequency is modified by the inhomo- 
geneous coupling. We remark that these oscillation fre- 
quencies are valid only for small amplitude oscillations 
around the potential minimum, the trap as a whole be- 
ing far from harmonic. 

For values of Wri/^i close to (3, one will face a situ- 
ation where the potential is almost flat in the crescent 
defined by z < 0, x = and 4y2 = (fig. [Jf,). The 
gravitational energy mgro at the holes is indeed exactly 
equal to the energy shift Fhfli induced by the rf at the 
bottom of the iso-i? surface, and ujy as expressed above 
goes to zero. However, this intermediate region is ex- 
tremely narrow for /3 ^ 1, with a width in the parameter 
Ljy[/Qi on the order of 3/(2/3) around /3, and in practice 
one deals only with one of the two extreme situations, 
a or c. Note that in any case the zero magnetic field 
region, responsible for Majorana spin flips in the centre 
of conventional quadrupole traps, plays a minor role in 
an rf-dressed quadrupole trap. The energy barrier is at 
least Fhujj-f, much larger than the thermal energy of the 
atomic cloud. 

In our experiment, ajj-f/fli is larger than (3 (typically 
ujr:i/2'K ~ 3 MHz and Q.i/2n < 260 kHz), and we are 
always in the case of a single minimum. Therefore, due 
to gravity, ultracold atoms remain trapped near the bot- 
tom of the ellipsoid and only marginally explore the non- 
coupling regions. This single- well situation was demon- 



strated experimentally in a loffe-Pritchard type trap [5| . 
The opposite situation was used in an atom chip exper- 
iment for producing a double well potential frj. In the 
experiment described in this paper, we bring into light 
the non negligible effect of the two holes for an ultra- 
cold atomic cloud, even in the situation where a single 
minimum is present. 



III. EXPERIMENTAL SEQUENCE 

The experimental setup has been described in detail 
elsewhere [l^ The atoms in the state F = 2, 

niF — +2 are first loaded into a static magnetic trap, in 
a Quadrupole and loffe Configuration (QUIC) see 
fig. m Quad and loffe coils. In this trap, the atoms are 
confined near the magnetic field minimum, position A of 
fig.m situated 7 mm away from the quadrupole center O 
in the direction of the loffe coil. Then, a 30 s radio fre- 
quency ramp is applied to evaporatively cool them down 
to just above the condensation threshold. We deliber- 
ately do not reach condensation and rather work with 
ultracold clouds, with a temperature of typically 4 /iK, 
in order to explore the leakage through the pierced part 
of the iso-B bubble as explained later. The evapora- 
tion rf source is then switched off and the trapping rf 
field, polarized along y, is turned on at a frequency lower 
than the resonance frequency in the trap centre, 1.3 MHz. 
Note that for this purpose, a non zero minimum magnetic 
trap is necessary. The rf frequency is then adiabatically 
ramped up to the desired value of uj^{ within typically 
150 ms, as detailed in This way, the atoms are al- 
ways following the upper dressed state, corresponding to 
niF = F = +2 in the dressed basis. At the end of this 
step, they are confined in a dressed rf trap relying on the 
QUIC static magnetic field. 

The QUIC trap stage is necessary to obtain an ultra- 
cold sample and load it efficiently into the dressed rf trap. 
The atoms are then transferred into a dressed quadrupole 
trap, which presents points of strictly zero rf coupling at 
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any value of Wrf . The current in the lofFe coil is therefore 
decreased in order to reach a quadrupole configuration 
within a time interval of about 450 ms. During this pro- 
cess, the minimum of the static magnetic field quickly 
goes to zero and is split into two zero field points sepa- 
rating along the x axis. One point is going back 7 mm 
away from A to the position O of the initial quadrupole 
centre whereas the other one is going from A to infinity 
in the opposite direction (see Ref. Fig. 2). During 
this transfer, the atomic cloud remains in the dressed 
trap, situated just below the static magnetic minimum, 
and is split into two clouds according to the deformation 
of the iso-magnetic surfaces. The first stage of the loffe 
current decreasing is controlled carefully in order that at 
least half of the atoms follow the right direction, toward 
the initial quadrupole centre O. At the end of this stage, 
about 4 X 10^ atoms at a higher temperature of 8 /iK are 
confined in an rf-dressed quadrupole trap. The QUIC-to- 
quadrupole transfer of the rf-dressed atoms is a critical 
step and is very sensitive to initial conditions. This pro- 
cedure is at the origin of a scattering in the values of 
atom number and temperature. 

The rf field is produced by a home made synthe- 
sizer [1^. The Rabi frequency fii corresponding to a 
given voltage amplitude is evaluated as follows. The rf 
antenna, normally dedicated to rf trapping, is used for 
evaporative cooling in the QUIC trap for this measure- 
ment. An evaporation ramp is applied down to a given 
final rf frequency. The threshold frequency below which 
all atoms are expelled from the magnetic trap is recorded. 
For an arbitrary weak rf power, this threshold frequency 
corresponds to the Zeeman splitting at the bottom of 
the QUIC magnetic trap. However, the recorded thresh- 
old frequency is larger, as the magnetic levels are de- 
formed due to dressing by the rf photons. We use the 
observed value as input in the calculation of the dressed 
potential, including gravity, and search for the Rabi fre- 
quency value making the potential flat, just unable to 
hold the atoms. We repeat this procedure for different 
voltage amplitudes at the synthesizer to calibrate the rf 
power. 



IV. EXPERIMENTAL RESULTS 

After the atoms are transferred to the rf-dressed 
quadrupole trap, they remain stored for a variable time 
T after which the magnetic field and rf coupling are si- 
multaneously switched off. An absorption image of the 
cloud is taken after a 7.5 ms ballistic expansion. From 
these measurements, the evolution of the atom number 
and the temperature with the storage time r is deduced. 
The temperature is related to the vertical size after ex- 
pansion of the atomic cloud. In principle, it also depends 
on the initial vertical size, which is not entirely negligible 
due to the curved shape of the trap. However, we neglect 
the initial size, which overestimates the real temperature 
by about 20%, as deduced from the comparison with a 




FIG. 4: Temperature T, measured as described in the text, 
as a function of time with f2i = 27r x 40 kHz and for different 
values of the RF frequency Wrf/Syr = 2 MHz (closed circles) 
and 3.1 MHz (open squares). Each point is an average over 2 
to 3 data points. The error bars take into account the point 
dispersion as well as an experimental uncertainty of 10%. T 
drops over the first trapping seconds, and then stabilizes to a 
value T = 1.35 /iK (resp. 3.68 indicated by the horizontal 
lines, deduced from the average of the 5 (resp. 4) last points. 



numerical simulation, see next section. 

We observe a rapid decreasing of the temperature dur- 
ing the first few seconds. The experiment was repeated 
at rJi/27r = 40 kHz for two different values of Wrf (see 
fig. 2]). For cjrf/27r ~ 2 MHz, we find a final tempera- 
ture Ti — 1.35 fxK whereas it does not get lower than 
T2 = 3.68 /.tK for u},i/2-K = 3.1 MHz. For the largest rf 
coupling value f2i/27r = 260 kHz, the temperature de- 
crease, can be observed over more than ten seconds (see 
the experimental points on fig. [5]), with a first drop in 
temperature followed by a transition to a slower but still 
decreasing evolution. 

These results, together with the non exponential de- 
crease of the atom number, see inset of fig. [5l reveal an 
evaporation process via Landau-Zener losses through the 
holes along the rf coil axis in the equatorial plane. The 
final temperature is related to the trap depth, calculated 
between the trap bottom and the position of the holes 
and depending on the rf frequency through 



a 



(9) 



The trap depth is Ua/kB — Qi — 26.7 /iK (resp. 
62 = 15.7 ^iK) for a 3.1 MHz (resp. 2 MHz) rf frequency 
and a Rabi frequency rii/27r = 40 kHz. The final tem- 
perature compared to the trap depth gives Qi/Ti = 7.3 
and Q2/T2 = 11.6. The overestimation of temperature 
makes these values a lower bound to the depth to tem- 
perature ratio. A ratio of 8 to 10 is considered typical 
for a spontaneous evaporation process in an harmonic 
trap 

In principle, the evaporation could also come from 
losses through the centre of the quadrupole, however the 
energy Fharo = FTioj-cf required to cross the central re- 
gion of the quadrupole trap is much higher than the en- 
ergy Uq necessary to explore the uncoupled region, mak- 
ing this process very unlikely. 
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FIG. 5: Temperature, deduced from the size after a 7.5 ms 
time-of- flight, as a function of the holding time in the rf- 
dressed quadrupole trap. Experimental data, closed dia- 
monds with error bars as in Fig. |4j are compared to the nu- 
merical simulation, thick line. Inset: Evolution of the atom 
number in the same conditions: experimental data, circles, 
are compared to the simulation, line and to a pure exponen- 
tial fit, dashed red line. For these experiments and the re- 
lated simulation, the parameters are uiii/2n = 3.1 MHz and 
^1 = 27r X 260 kHz. 



Finally, the possible non-adiabatic Landau-Zcncr 
losses in the avoided crossing region at the bottom of the 
ellipsoid, where the atoms mostly sit, can be neglected. 
Indeed, the probability of spin-flipping in this zone for a 
given velocity v reads [ist 



Plz = 1- 



1 - cxp ^ 

^ ' 2av 



n 2F 



(10) 



If one averages Eq. (flOl) over the whole velocity distri- 
bution for an rf coupling ili = 30 kHz, smaller than was 
ever used in the experiment, and a temperature in the 
dressed trap T = 10 /xK, one finds a life time tlz due to 
Landau-Zener losses as high as 100 s, which is even larger 
than the life time due to collisions with the background 
gas. Landau-Zener losses at the trap bottom are then 
obviously not responsible for the losses and temperature 
decrease. 



V. NUMERICAL SIMULATION 

In order to confirm our interpretation, we performed 
a numerical simulation of the cloud classical dynamics 
with the parameters of the experiment, for a coupling 
f]i/27r = 260 kHz, a rf frequency ui^{/2tt = 3.1 MHz 
and a trapping time ranging between and 10 s. The 
simulation calculates the evolution of a cloud of atoms 
initially trapped on the iso-B surface, at temperature 
Ti, above condensation threshold. Collisions between 
trapped atoms are taken into account through a method 
developed by G.A. Bird [l^, while collisions with the 
background gas are described through a lifetime of 60 s. 
A possible linear heating rate H is taken into account 
in the simulation by increasing each velocity component 



FIG. 6: Estimated phase space density for the data series 
of Fig. [5l from the experimental data, diamonds, and the 
simulation, line. The parameters are ujri/2n = 3.1 MHz and 
i}i — 2n X 260 kHz. The error bars take into account data 
dispersion and an experimental uncertainty of 40%, resulting 
from uncertainties in atom number and temperature. 



Vi by 2viHdt/T after a time step dt. The Landau-Zener 
losses are more difficult to model. The use of a random 
number to decide, with a local probability distribution, if 
an atom is lost or not would be time-consuming. Instead, 
the holes are modeled by a deterministic loss term: the 
atoms are considered lost if the Landau-Zener probability 
at this point Plz{y,v) exceeds a critical value Pc- Due 
to the exponential dependence of the Landau-Zener non 
adiabatic coupling on the rf Rabi frequency, this occurs 
essentially in a region very close to the two holes. 

Fig.[5]gives a plot of the temperature for different hold- 
ing times in the dressed trap. The temperature is de- 
duced in the same conditions as in the experiment, from 
the cloud size after a 7.5 ms time-of- flight, and the two 
values are thus directly comparable. The simulation flts 
at best the experimental data for the choice Pc ~ 0.9 and 
H = 300 nK/s. The Landau-Zener probability Pc mainly 
determines the initial rate of temperature decrease. Once 
this value is fixed, the heating rate H imposes the final 
saturation temperature, resulting from the equilibrium 
between evaporation and residual heating. The simula- 
tion reproduces well the features of spontaneous evapora- 
tion, namely temperature reduction and non exponential 
atomic losses, see dotted line in the inset of fig. [5l These 
features completely disappear if the holes are suppressed. 

The peak phase space density in the rf-dressed trap is 
computed from the simulated data. To be able to com- 
pare this figure with the experiment, we used in both 
cases the estimation Lp ~ N{hOj/kBT)^, where T is cal- 
culated from the vertical width after time-of-flight as ex- 
plained above. Here, (D is the geometric mean of the os- 
cillation frequencies of the rf-dressed trap, Eq.([6]l8|). This 
estimate only gives the order of magnitude for the phase 
space density, as the formula is correct only for a har- 
monic trap and the trap is far from being harmonic. The 
results are displayed on Fig. [SI for the same parameters 
as used in Fig. [51 For these parameters, ui/2tt = 60.6 Hz. 
The simulation shows a flrst increase of about 50 % 
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within the first 3 seconds, followed by a regular decrease 
due to residual heating. The experimental data give val- 
ues for if of the same order, but are too scattered for 
demonstrating unambiguously a phase space density in- 
crease. This is due to the sensitivity of the estimate for 
ip to the temperature value T, which is determined with 
a large uncertainty. 

VI. CONCLUSION 

In this paper, we demonstrate trapping of ultracold 
atoms in a quadrupole magnetic configuration. Atoms 
are prevented to escape through the centre by rf-dressing 
of the atomic state, which results in a confining adia- 
batic potential. Non exponential losses associated to a 
temperature decrease are observed in this trap, giving 
evidence for spontaneous evaporation. These losses are 
attributed to the presence of zero rf-coupling zones in the 
trap, due to rf polarization. Experimental data are inter- 
preted with the help of a numerical simulation, well re- 
producing the main features. The evaporation results in 
a drop in temperature, which saturates due to a residual 
heating in the trap. This technical heating clearly pre- 
vents an increase in phase space density. Cooling in the 



trap may be improved in different ways. The observed 
spontaneous evaporation could be forced by reducing the 
rf frequency ojif with time, thus limiting the trap depth. 
Still, it seems that a necessary approach would be to 
use a second independent rf field to spin flip selectively 
the more energetic atoms, as in conventional evaporative 
cooling [2^, l2l!|. Ultimately, the quadrupole geometry 
with a symmetry axis oriented in the vertical direction is 
ideal for the realization of a isotropic quasi-2D degener- 
ate gas jl^l or an atomic ring Q. 
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